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Abstract 

We investigate the full counting statistics of a voltage-driven normal metal(N)-superconductor(S) contact. In the 
low-bias regime below the superconducting gap, the NS contact can be mapped onto a purely normal contact, albeit 
with doubled voltage and counting helds. Hence in this regime the transport characteristics can be obtained by the 
corresponding substitution of the normal metal results. The elementary processes are single Andreev transfers and 
electron- and hole-like Andreev transfers. Considering Lorentzian voltage pulses we hnd an optimal quantization for 
half-integer Levitons. 
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1. Introduction 

Quantum shot noise and full counting statistics (FCS) 
have emerged as central tools of quantum transport in 
the last two decades. The main driving force is the dra¬ 
matic difference in properties of fluctuations of the cur¬ 
rent of classical particles versus quantum particles be¬ 
having sometimes in a wave-like fashion |[I]. Classi¬ 
cal particles in a tunneling setup lead to fluctuations in 
the current described by Schottky’s formula for Poisson 
noise S - el, where S is the noise power of current 
fluctuations, I is the average current and e is the elec¬ 
tron charge ||2|. Considering the wave-like nature of 
electrons, which is encountered in nanostructured con¬ 
ductors at low temperatures, in combination with the 
quantum statistical fermionic Pauli principle leads to a 
suppression of the shot noise by a so called Fano factor 
F - Yin A(1 “ Tn)l Y Tn, where r„ are the transmis¬ 
sion probabilities of the electron waves in channels n 
EH). The suppression has been experimentally veri¬ 
fied in quantum point contacts 00 and other coherent 
conductors like diffusive wires with the characteristic 
Fano factor F = 1/3 El[8]I). 

Another leap forward in the understanding of quan¬ 
tum transport was to go beyond the average current 
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and the noise by considering the full counting statis¬ 
tics (FCS) of the transferred charge, which comprises 
all probabilities P(V) to transfer N charges. Equiv¬ 
alently, one considers the cumulant generating func¬ 
tion (CGF) SOtf) - ln[{e'^^)]. The remarkable result 
for a quantum contact at low temperature is SC\^) - 
{leVtalh) Yn lit [l + T„(e'^ - l)j which describes a bi¬ 
nomial distribution for each channel ifTOl . Note that 
in general a decomposition of the CGF into binomials 
or multinomials allows to identify the elementary pro¬ 
cesses and their probabilities. 

In superconductors the electrons are correlated in a 
single macroscopic wave function, which describes a 
condensate of so-called Cooper pairs consisting of two 
bound electrons. The state is stabilized by a finite bind¬ 
ing energy A, which needs to be payed twice to break up 
a Cooper pair into two independent electrons. In quan¬ 
tum transport this is manifest in an energy gap A be¬ 
low which the differential conductance of a junction be¬ 
tween a normal metal and a superconductor due to sin¬ 
gle electrons vanishes. However, the electron transport 
is still possible by an intriguing process called Andreev 
reflection in which a Cooper pair is transferred into a su¬ 
perconductor while the hole-like quasiparticle is left be¬ 
hind iini- Hence, in this process two charges are trans¬ 
ferred which is therefore possible also at subgap ener¬ 
gies. This process occurs with the probability of An- 
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dreev reflection= T^l(2- T„f- ||T2l[l3l. It is very 
interesting to note that the FCS for Andreev reflection 
takes a very similar form as in the normal case, namely 
^SaCt) = (2eVtolh) Zn In [l + - 1)] M- There¬ 

fore, the statistics is also binomial, but with the im¬ 
portant difference that in each process 2 charges are 
transferred. This follows from the 7r-periodicity due 
to a doubling of the counting field x in the factor 
_ I 'pjjg doubling of the effective charge trans¬ 
ported is manifest in the ratio between noise and av¬ 
erage current SII - eF^, where the Fano factor is now 
Fa = ^YjnRnil -R^)IZnRt In particular. Fa = 2/3 
in a diffusive normal-metal - superconductor junction 

nsiiihiini. 

Time-dependent voltage drives can be used to probe 
the dynamics of electrons in transport. One interesting 
aspect is that with a signal with a finite frequency to 
one has a tool to access the internal time-scale of the 
manybody state given by eVjh at low enough temper¬ 
atures. This shows up, for example, in the noise of a 
quantum contact driven by harmonic voltage. The noise 
is a piecewise linear function of the dc voltage bias with 
slopes which depend on the amplitude of the ac voltage 
component. The kinks in noise occur at dc bias voltages 
eV - nhto matching an integer multiple of the drive 
frequency US). In the normal case, the noise in the 
presence of the drive is always larger than or equal to 
the dc noise level. The scattering theory of the excess 
photon-assisted noise has been put forward by Pedersen 
and Biittiker ifT^ . The further advancement was inter¬ 
pretation of the noise and current cross-correlations in 
terms of excited electron-hole pairs that was given by 
Rychkov, Polianski, and Biittiker 1201 for an ac drive of 
low amplitude, eVac Hoo, where at most one electron- 
hole pair can be created per voltage cycle. Remarkably, 
this picture of electron-hole pairs created by the drive 
persists even at large amplitudes and to all orders in 
charge transfer statistics 112111221 . Photon-assisted noise 
has been observed experimentally in normal coherent 
conductors EH El and in diffusive normal metal - su¬ 
perconductor junctions 1251 . More recently, quantum 
noise oscillations have been observed in a driven tun¬ 
nel junction l26l . The noise spectral density for dc- and 
ac-bias voltages for normal metal-superconducting con¬ 
tacts has been discussed in l34l . 

An extremely intriguing possibility is the ability to 
control the electron dynamics by shaping the voltage 
pulses. In particular, it was shown that Lorentzian volt¬ 
age pulses with a quantization condition e f dtV(t) - nh 
result in the soliton-like electronic excitations which 
minimize the noise level to the one of an equivalent dc 


voltage, S = e/F lE2l|28]. These so-called levitons 
are hence a collective single-electron excitations local¬ 
ized in space and time, which offer interesting perspec¬ 
tives as carriers of quantum information Il29l [30l . To 
access the full counting statistics in the presence of a 
time-dependent drive a non-equilibrium quantum held 
theoretical approach to FCS was developed by Nazarov 
and one of the authors ED. This allowed to perform the 
analysis of the FCS in terms of elementary events for an 
arbitrary time-dependent voltage V{t) ED ED- The re¬ 
sults is that one has to distinguish two types of events: 
Single electron transfers, which occur with a frequency 
of the average voltage and have the standard binomial 
statistics = (2eydcfo//j) In [l + - 1)] 

and electron-hole pairs obeying a trinomial statistics 
Sehix) = Y,nk Mk In [1 H- 2T„R„pk(cos(x) - 1)]. The 
probabilities ph are interpreted as probabilities of 
electron-hole pair creations and depend in a character¬ 
istic way on the driving voltage which is assumed to be 
periodic with frequency to. The number of attempts for 
the pairs to traverse the contact is M oc 2htotQ/h. This 
opens a route towards dynamic control of elementary 
excitations using suitably tailored voltage pulses ED- 

In this article, we consider the FCS of an Andreev 
contact driven by a time-dependent voltage. Using an 
exact mapping of an NS contact onto an effective nor¬ 
mal contact E3^ we identify the elementary Andreev 
events and characterize the two types of processes. Sin¬ 
gle Andreev-pair transfers have binomial statistics and 
are determined by the average dc voltage Vdc. The time- 
dependent drive manifests itself in correlated electron- 
hole pairs which are transferred coherently. The respec¬ 
tive probabilities are found from the normal ones by the 
mapping p^[V(t)] «-> p^{2V{t)\. Indeed, by consid¬ 
ering as an example the Lorentzian voltage pulses we 
hnd a maximal noise suppression for half-integer pulses 
e f V(t)dt - nhl2 with integer n. Furthermore, increas¬ 
ing the voltage level in an ac-driven contact above the 
gap, we find a transition to minima at integer quantized 
voltages eVijc = nhto. 

The article is organized as follows. In Sec.[^ we in¬ 
troduce the extended Keldysh Greens function theory of 
quantum transport applied to a time-dependent voltage 
drive. In Sec.|^ we obtain the mapping of an NS contact 
to an NN contact and analyze the resulting FCS in terms 
of elementary events. Finally, in Sec. I^we discuss some 
examples of a voltage drive and consider the transition 
from Andreev to normal transport for large biases. 
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2. Keldysh formulation of Andreev contacts 


The Keldysh Greens function formalism is a very 
powerful method suitable for quantum nonequilibrium 
problems. We formally introduce the standard closed 
time-path and define Greens functions on the contour 
G(t,t') = mapped onto Keldysh space. 

Treating the time variables on the upper and the lower 
branches of the contour as independent, one can define 
a matrix Greens function 




( 1 ) 


In the quasiclassical approximation for a free Fermi gas 
at equilibrium, the Greens function reads 


G{e) - -inNo 


I l-2/(e) 
^(l-/(e)) 


-2/(e) 'l 
2/(e)-lj 


( 2 ) 


The prefactor Nq containing the density of states at the 
Fermi level is usually removed by proper normalization, 
so that the Greens function obeys the normalization con¬ 
dition = 1. Adding superconductivity results in the 
replacement -> (ij/e{t),il/h{t)Y = 
and hence an extension of the Keldysh matrix space by 
an additional electron-hole degree of freedom (Nambu 
space). 

For the present purpose a representation is chosen in 
which the Keldysh matrices are blocks in the Nambu 
space. Flence, the Keldysh Greens function of the nor¬ 
mal lead is given by: 

_y. m 

where check('') denotes matrices in Nambu(“) ® 
Keldysh(') space. Here, we introduce the (Keldysh- 
rotated) Greens functions 




1 2UhW 
0 -1 


_ /I 2UUiU\ 


(4) 


where h - I - 2f, f{e) - [exp(/3e) + 1 ]“' is the 
Fermi distribution function, and U is related to the time- 
dependent drive y(f), 

U(t', t") = e-'/o evmih ( 5 ) 


To access the FCS the counting field is incorporated 
into the Green’s function as 

Gn(x) = (6) 
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where f/f = T 3 ® fi. This gives 


Gn(x)^ 




( 7 ) 


The Greens function of the superconducting lead at 
low temperatures and drive energies well below the gap 
(ksTe, \eV(t)\ « A) is given by 


Gs -T 2 ® I - 


0 
i i 



( 8 ) 


Cumulant generating function is given by ED: 

^SOr) = ^Trln(l + ^({G^O^),G^( 0 )) - 2 )), (9) 


where Tr stands for the trace in Nambu, Keldysh, and 
time (energy) indices, and also implies a summation 
over transport channels T„. [y'-independent constant 
which ensures S(x = 0 ) = 0 is omitted for brevity.] 
Note that in deriving Eq. (|^ it was assumed that the 
dwell time in the scattering region is very short, so that 
the energy dependence of the scattering amplitudes can 
be neglected. 


3. Full counting statistics analysis 

At low energies, when Gs is given by Eq. ([ 8 ]), the 
CGE reduces to the normal-state circuit with the elec¬ 
tron and the hole Green’s functions: 

= ^Trln^l + ?i[{G%(x),GU-x)} - 2)j . (10) 

Note that now G^* are the normal-state Green’s func¬ 
tions in the Keldysh space, and the transmission proba¬ 
bilities r„ are replaced by Andreev reflection probabili¬ 
ties, 7?]^ = f{2-T„)^. After carrying out a gauge trans¬ 

formation, it is possible to ascribe the counting field and 
the drive to one ’lead’ only, and we obtain 

^SCr) = ^Trlnjl + ^({G^(2y(f),2;r),G“) - 2)j , 

( 11 ) 

where G^° = G^(F = 0,x - 0). Thus, at the sub¬ 
gap energies, the system can be mapped to a normal- 
state circuit with a doubled voltage drive and a doubled 
counting field, and with transmission probabilities T„ 
replaced by R^. This mapping is shown in Eig.[^ 

CGE in Eq. o can further be brought in the form: 

*S0r)=Trln[l+(l-/)/<(c2'l^-l) 

+ /(l-/)<(e'"'^-l)], ( 12 ) 
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Figure 1: An Andreev contact (a) with transmission probabilities {r„} 
between a normal and a superconducting metal is mapped onto a con¬ 
tact (b) between two normal metals (electron and hole space), where 
the transmission probabilities are replaced by the Andreev reflection 
probabilities (c) By a further gauge transformation the counting 
field and the (time-dependent) voltage are applied on one side only. 


where / = accounts for an effective dou¬ 

bling of the drive voltage. At zero temperature the 
matrix operators /, / have the additional property that 
- f and = /, which allows us to decompose the 
PCS into the single electron processes and electron-hole 
pairs as mentioned in the introduction. Here in the An¬ 
dreev case they take a slightly different form 

•Si = ^^^^ln[l+f;;J(e'2^-l)], (13) 

n 

Seh = M In [l + 2Ri(l - R^)p,(cosdx) " O] ■ 

nk 

(14) 

CGF in Eq. ( |T3| ) accounts for the Andreev reflection 
of the excess electrons due to dc voltage applied, cf. 
Fig. I^a). The charges are transferred in pairs and the 
statistics is binomial in each transport channel with the 
transmission probabilities given by the Andreev reflec¬ 
tion coefficients R^. The rate leVadh with which the 
excess electrons impinge on the contact is the same as 



Figure 2: Elementary transport processes in a driven NS contact: (a) 
Andreev reflection of an excess electron due to dc voltage applied, see 
Eq. (b) Andreev reflection of an electron from the electron-hole 
pair accompanied by the normal reflection of the hole from the pair, 
see Eq. O- A reverse process is also possible in which the hole ex¬ 
hibits Andreev reflection and the electron exhibits normal reflection. 
Electrons and holes are defined relative to the energy of the super¬ 
conducting condensate (6 = 0). Time-dependent drive is indicated by 
shading. 


in the normal case. CGF in Eq. Cl accounts for the 
charge transfer statistics due to electron-hole pairs cre¬ 
ated by the ac drive. In general, S^h consists of two 
types of electron-hole processes with different proba¬ 
bilities pk of the electron-hole pair creations and dif¬ 
ferent numbers of attempts for particles to traverse the 
junction, M - (hciJtolh)(l - v) and M - (/iwfo//z)v. 
Here, v = leV^^lhu - \2eVii^lfui)\ is a fractional part 
of 2eVdclfiCL> (L-^J denotes the integer part of x). Charge 
transfers described by Eq. d are bidirectional pro¬ 
cesses in which an electron from the pair is Andreev 
reflected at the contact while the hole exhibits a nor¬ 
mal reflection, or vice versa. This is schematically de¬ 
picted in Fig. I^b). As a result, pairs of charge quanta 
are transferred in either direction with the probabilities 
R^il - Rn)Pk in each transport channel. 


4. Andreev levitons in an NS junction 


As a first example we compute the photon-assisted 
noise when the junction is driven by periodic Lorentzian 
voltage pulses of the form; 


vu,m = ^ T 

71 


Ttl 


it-kTd + rl 
Vo sinh(27rT/,/r) 
cosh(27TTilT) - cos(27rf/7’) 


(15) 


Here, T - 27r/w is the period of the drive, and the am¬ 
plitude Vo is chosen in such a way as to represent the 
average voltage per period, {l/T) - Vo- 

Lorentzian voltage pulses are shown in Fig. for dif¬ 
ferent widths TijT - 0.05, 0.1, 0.2. 
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Figure 3: Lorentzian voltage pulses for different widths: t/,/ 7 = 0.05 
(dash-dotted), 0.1 (solid), and 0.2 (dashed line). 


The excess photon-assisted noise (the total noise mi¬ 
nus dc noise level) is given by 


c ” 

ac _ \ ^ 

J-L 


2eVo 

1 |2 

2eVn 

hcj ' " 

\an\ - 

hai 


(16) 


where 


5o = y<(l--R^). (17) 


Coefficients a„ are related to the doubled ac part of the 
drive voltage, 2yac(f) = 2(yLor(0 - yo): 




1 

T 







(18) 


where (p{t) - (elh) j^2V^c(t')dt'. These coefhcients 
read 


„y (-l)V"^r(g + « + ^) 

^ r{k + i)r(^ -k + i)r(n + k + i) 


(19) 


for n > 0, and 


i„i y (-i)^^”y2cr(g + ^) 

^ ^(k + l)r(^ -\n\-k+ l)r(|n| +k+l) 

( 20 ) 

for n < 0. Here, q - 2eVo/hu and y - . 

The photon-assisted noise is shown in Fig. |^a). As 
the width is increased, the pulses overlap more 
strongly and yLor(f) approaches the constant voltage 
Vq. This results in the overall suppression of the ex¬ 
cess photon-assisted noise. In addition, the excess noise 
is fully suppressed at half-integer values of eVo/hu. 
This corresponds to the half-integer Lorentzian pulses 



eV^/hw 


Figure 4: (a) Excess photon-assisted noise due to a Lorentzian voltage 
drive shown in Fig. and (b) the corresponding probabilities of the 
electron-hole pair creations as a function of the drive amplitude. 


e f VLoAOdt = nh/2, in contrast to the normal junc¬ 
tions in which the noise suppression occurs for integer 
pulses. The photon-assisted noise can also be expressed 
in terms of elementary events of the electron-hole pair 
creations EHES: 

S,JSo = 2(l-v)'^Pi. (21) 

k 

Here, v = 2eVolhu - \ 2eVQ jfiLol is the fractional part 
of 2eVo/ha) and pi^ are the probabilities of the pair cre¬ 
ations given by pk - sm^(akl2) where are the 
eigenvalues of the operator hh. The electron-hole cre¬ 
ation probabilities as a function of the amplitude of the 
Lorentzian pulses are shown in Fig. Qb). We find that 
in the problem at hand, there is only one electron-hole 
pair created per period with probability Pi s p. The pair 
creation probability increases as eVolhoj approaches the 
half-integer values. However, the photon-assisted noise 
is nevertheless zero at these points because the effective 
rate of attempts ha>(l - v)/h vanishes. 

5. Large excitation noise 

In what follows we allow for the dc bias and the ac 
drive amplitudes to be comparable to A. For simplicity, 
we still assume low-temperature limit, Tg - 0. 
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When the drive amplitudes are comparable to A, the 
NS junction can no longer be mapped to the normal one. 
However, we can proceed with the numerical calcula¬ 
tion of the cumulant generating function in Eq. Q. The 
Green’s functions are given by 


Giv(O) = 


T3 

0 



/ Gr {Gr - GA)h{e)\ 

i 0 Ga j- 


( 22 ) 

(23) 


Here, 



Gr,a - 


+ 1 


V(e + i0)2 - A2 


e ± iO 


-{e±iO) 


(24) 

(25) 


and hi - UhU\ /i2 = U''hU. 

Next, we note that for the periodic time-dependent 
drive with the period T - Injto, the operators hi 2 cou¬ 
ple only the energies that differ by an integer multiple 
of ha>. Therefore we can use a matrix representation 
in energy indices, {hi)nm{G} - hi{e + nha>, e + mhcj) 
{-ha>l2 < e < fi.aj/2). This provides a matrix structure 
in energy in Gat and Gg. The trace operation in Eq. 
now amounts to a matrix diagonalization in Keldysh, 
Nambu, and energy indices and integration over e. 

We calculate the cumulant generating function for a 
diffusive NS junction that has a distribution of transmis¬ 
sion eigenvalues given by 


P(T) = 


Gn 1 

Gq 2^ Vi - t’ 


(26) 


where Gfj is the normal-state conductance and Gq - 
2e^/h and it is assumed that the Thouless energy is 
much larger than all other relevant energy scales. The 
cumulant generating function in Eq. 0 reduces to 

>S(;r) = - ^Tr arcsin^ |i ^2 - {G^CrXGslj ■ (27) 

Eor the average current and the current noise power we 
obtain 


I = G, 


de d arcsin^( gjA(e,x)l2) 1 


J-hu/2 4e ^ 


S =-Ga 


ddx) lo 
de arcsin^C gjA(e,x)l2) 


(28) 


de y 

J-ruj/2 4 


d(ix)^ 


(29) 


where A are the eigenvalues of 2 - {GA'Ct')> G, 


SI- 



Figure 5: Excess photon-assisted noise for the Lorentzian pulses with 
tlIT = 0.1. The noise is shown for different drive frequencies: 
/ku/A = 0.25, 0.20, 0.15, 0.10, 0.02 (top to bottom). 


Photon-assisted noise 5ac - S - Sic for an NS junc¬ 
tion driven by periodic Lorentzian pulses is shown in 
Eig. I^for different drive frequencies /Uu/A = 0.25, 0.2, 
0.15, 0.1, 0.02 (top to bottom). The noise is normal¬ 
ized to 5 0 in Eq. which in the case of a diffusive 
junction reads So = G^hajp. 

The half-integer Lorentzian pulses with eVolfdo - 
n/2 at energies much smaller than the superconduct¬ 
ing gap create Andreev levitons which are the minimal 
excitation states of the NS system. As the energy be¬ 
comes comparable to the gap, the system can no longer 
be mapped to an effective normal junction. Therefore, 
the half-integer Lorentzian pulses are no longer optimal 
due to contributions of normal and Andreev electron 
transport above the gap. As a result, the excess photon 
assisted noise starts to increase, see Eig. However, 
at energies much larger than the gap, the junction is in 
the normal state and the integer Lorentizan pulses create 
minimal excitation states in the normal junction. 

The quantum oscillations of the photon assisted noise 
as a function of the dc voltage have been observed re¬ 
cently in the normal-state tunnel junction driven by har¬ 
monic time-dependent voltage y(f) = yo[l + cos(a;i)] 
ES\ . Harmonic drive in general creates additional 
electron-hole pairs which results in the non-zero excess 
noise with minima at integer values eVolhcj - n. The 
same is true in the NS junction at energies much lower 
than the gap, except that the excess noise minima appear 
at half-integer values eVolhoj - n/2, see Eig. At in¬ 
termediate drive frequencies that are comparable to the 
gap, the noise in the system is determined by a density 
of states which is affected both by the drive and by the 
superconducting proximity effect. In addition, both nor¬ 
mal and Andreev processes contribute to the transport. 
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Figure 6: Excess photon-assisted noise for the harmonic time- 
dependent voltage V(t) = Vo[l + cos(6j?)] and amplitudes of the 
drive that are (a) smaller and (b) larger than the gap. At subgap en¬ 
ergies the excess noise is minimal at half-integer drive amplitudes 
eValhcj = njl, while above the gap it is minimal at integer drive 
amplitudes eVolhcj = n. 

As a result, the total noise can even be suppressed be¬ 
low the effective dc noise level (the excess noise can be 
negative), the situation which otherwise cannot occur in 
the normal junction. 

6. Conclusion 

We have analyzed the transport properties of a driven 
quantum point contact between a normal metal and a 
superconductor. Using an extended Keldysh Greens 
function method we could determine the full counting 
statistics and identify the elementary transport processes 
for an arbitrary voltage drive in the subgap regime. At 
voltage amplitudes and frequencies well below the su¬ 
perconducting gap A the NS-contact can be mapped 
onto a contact with two normal leads. The transmis¬ 
sion is determined by the Andreev reflection probability 
and the probabilities of elementary Andreev processes 
are determined by the ones known for normal transport 
with an effective charge 2e. In that spirit we have dis¬ 


cussed Lorentzian voltage pulse which lead to Andreev- 
Levitons, which are now pure two-charge excitations, 
for half-integer quantized amplitudes e J dtV — nhjl. 
Finally we have discussed the transition from half¬ 
integer steps eVolho) - n/2 in the ac-noise to sim¬ 
ply quantized steps for voltages and frequencies much 
larger than A. 

In future it will be interesting to investigate open 
questions, e.g., for which parameters the noise is min¬ 
imized and the nature of the elementary events at in¬ 
termediate frequencies where both Andreev and normal 
reflection processes coexist, as well as the effects of de¬ 
phasing related to the finite Thouless energy. 

We acknowledge financial support by DFG through 
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